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Abstract 

In this article, we present the scalar-diquark-scalar-diquark-antiquark type and scalar- 
diquark-axialvector-diquark-antiquark type pentaquark configurations in the diquark model, 
and study the masses and pole residues of the = \ hidden-charm pentaquark states 
in details with the QCD sum rules by extending our previous work on the = | and 
hidden-charm pentaquark states. We calculate the contributions of the vacuum conden¬ 
sates up to dimension-10 in the operator product expansion by constructing both the scalar- 
diquark-scalar-diquark-antiquark type and scalar-diquark-axialvector-diquark-antiquark type 
interpolating currents. The present predictions of the masses can be confronted to the LHCb 
experimental data in the future. 
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1 Introduction 

Recently, the LHCb collaboration observed two pentaquark candidates Pc (4380) and Pc (4450) in 
the J/V'P mass spectrum in the A° —>■ J/ipK~p decays with the significances of more than 9tT 
[I]. The measured masses and widths are Mp^(43gQ) = 4380 ± 8 ± 29MeV, Mp^(445o) = 4449.8 ± 
1.7±2.5MeV, rp^(438o) = 205± 18±86MeV and rp^(445o) = 39±5± 19MeV, respectively. The 
Pc(4380) and Pc(4450) have the preferred spin-parity = | and I"*", respectively. The decays 
Pc (4380) J/ij^p take place through relative S-wave while the decays Pc (4450) J/ijjp take place 
through relative P-wave, the decays Pc(4450) —>■ J/ij^p are suppressed in the phase space, so the 
Pc(4450) has smaller width. There have been several attempted assignments, such as the EcP*, 
S*P>*, XciP, J/V'-^(1440), J/V'./V(1520) molecule-like pentaquark states [2] (or not the molecular 
pentaquark states m, the diquark-diquark-antiquark type pentaquark states [HE], the diquark- 
triquark type pentaquark states |6], re-scattering effects |7], etc. We can test their resonant nature 
by using photoproduction off a proton target [5]. 

In Ref. |5| , we construct the scalar-diquark-axialvector-diquark-antiquark type interpolating cur¬ 
rents, calculate the contributions of the vacuum condensates up to dimension-10 in the operator 
product expansion, and extend the energy scale formula suggested in our previous works [9] to 
study the masses and pole residues of the = | and | hidden-charm pentaquark states with 
the QCD sum rules, and assign the Pc(4380) and Pc(4450) to be the | and I"*" pentaquark states, 
respectively. In this article, we extend our previous work to study the diquark-diquark- 

antiquark type hidden charm pentaquark state by calculating the contributions of the vacuum 
condensates up to dimension-IO, and try to obtain the lowest masses based on the QCD sum rules. 

The article is arranged as follows: we choose the optimal pentaquark configurations in Sect.2; 
in Sect.3, we derive the QCD sum rules for the masses and pole residues of the ^ pentaquark 
states; in Sect.4, we present the numerical results; and Sect.5 is reserved for our summary and 
discussions. 
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2 Pentaquark configurations in the diquark model 

The diquarks qJCTq'f. have five structures in Dirac spinor space, where CT = Cqs, C, 

and Ccr^^ for the scalar, pseudoscalar, vector, axialvector and tensor diquarks, respectively, 
and the j and k are color indexes. The matrices and Ca^i, are symmetric, the matrices 
Cys, C and Cy^ys are antisymmetric. The attractive interactions of one-gluon exchange favor 
formation of the diquarks in color antitriplet 3c, flavor antitriplet 3/ and spin singlet Ig [10) . while 
the favored configurations are the scalar diquark states (e'-^^Q'J'Cysg^) and axialvector diquark 
states [TT1|T2]. The calculations based on the QCD sum rules indicate that the 

heavy-light scalar and axialvector diquark states have almost degenerate masses while the 
masses of the light axialvector diquark states lie about (150 — 200) MeV above that of the light 
scalar diquark states [T^], if they have the same quark constituents. In this article, we take the 
diquark states as basic constituents, and choose the scalar-diquark-scalar-diquark-antiquark type 
and scalar-diquark-axialvector-diquark-antiquark type pentaquark configurations. 

Now we illustrate how to construct the pentaquark states in the diquark model according to 
the spin-parity J^, 
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where the 1“ denotes the contribution of the additional P-wave to the spin-parity, the subscripts 
ud, uc, c and uudcc denote the quark constituents. The quark and antiquark have opposite parity, 
we usually take it for granted that the quarks have positive parity while the antiquarks have 
negative parity, so the c-quark has J+ = 4 . 

The overlined states | _ and I"*" . - are assigned to be the pentaquark states Pc(4380) and 

Pc (4450), respectively [S]. In previous work we choose the scalar-diquark-axialvector-diquark- 
antiquark type currents Jfj,{x) and 

Jfj.{x) = e"'‘°‘e"^'"e‘"""^uJ{x)Cj5dk{x)u^{x)Cjfj.Cn{x)Ccl{x), (5) 

Jfiuix) = -^s''‘°-e'^^''e^"^'^uJ{x)Cj5dkix)[u'^{x)C-ff^Cn{x)j„Cc'^{x)+u'^{x)CjuCn{x)j^,Cc'^{x)] , 

( 6 ) 

to interpolate the | and | pentaquark states, respectively, where the i, j, fc, • • • are color 
indices, the C is the charge conjugation matrix. 

The underlined states are supposed to be the lowest pentaquark states, while their P- 

wave partners are supposed to be the lowest pentaquark states with the positive parity. In 

this article, we choose both the scalar-diquark-scalar-diquark-antiquark type and scalar-diquark- 
axialvector-diquark-antiquark type currents 
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to study the lowest pentaquark states with in a consistent way, where the subscripts Jl 

and Jh denote the spins of the light and heavy diquarks, respectively. 


3 QCD sum rules for the pentaquark states 

In the following, we write down the two-point correlation functions (p) QCD sum 

rules, 

^ J d^xe^P-^0\T \0}. (9) 


The currents have positive parity, and couple potentially to the i"*" hidden-charm pen¬ 
taquark states , 

( 0 |J,.,h( 0 )|P+,„(p)) = A+^.„17+(p,s), ( 10 ) 

Kl3h are the pole residues, the spinors U^{p, s) satisfy the Dirac equations = 

0. On the other hand, the currents Jji^jn (0) also couple potentially to the ^ hidden-charm pen¬ 
taquark states as multiplying ^75 to the currents Jj^jh {x) changes their parity [I^[MlITKlllB] . 


muHmp,~,jp)) = KuHn5u-{p,s), (11) 

the spinors U^{p,s) (pole residues have analogous properties. 

We insert a complete set of intermediate pentaquark states with the same quantum numbers as 
the current operators JjLjnix), and into the correlation functions ^jLinip) to obtain 

the hadronic representation [niiiH]. After isolating the pole terms of the lowest states of the 
hidden-charm pentaquark states, we obtain the following results: 
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where the are the masses of the lowest pentaquark states with the parity ± respectively. 

We have to include the negative parity pentaquark states as _ according to the 

special quark configurations, see Eqs.(l-4). 

Now we obtain the hadronic spectral densities through the dispersion relation. 


^nilljz.ju (g) 

TT 


At... 


'jhjH “ + I ' ''3L3H V“ "'"3L3H, 

^^3L3H,+ ^tL3H (s “ ^jL3H, + ) ~ ^0L3H-^3L3H 

= ^^P]L3 His)+ 


(13) 


then we introduce the weight function exp (—^) to obtain the QCD sum rules at the hadron side 

/ Wsp]l3h (s) - P°L3H (s)] exp ( 

J 4m2 
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where the sq are the continuum threshold parameters and the are the Borel parameters. We sep¬ 
arate the contributions of the negative-parity (positive-parity) pentaquark states from the positive- 
parity (negative-parity) pentaquark states explicitly. 
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In the following we briefly outline the operator product expansion for the correlation functions 
^jLinip) in- perturbative QCD. Firstly, we contract the u, d and c quark fields in the correlation 
functions (p) with Wick theorem, and obtain the results: 

n„.(p) ^ 

{Tr [j,Dkk'ixh,CUj'^,ix)C] Tr [j,Cr.r.'ixh,CUZm'ix)C] 

-Tr [j5Dkk'ix)j5CU^j,(x)Cj5Cun'ix)j5CUj^,(x)C] 75 CCJ,^{-x)Cj5} , (16) 
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where the Uij{x), Cij{x) are the full u, d and c quark propagators respectively {Sij{x) = 
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and t” = y, the A” is the Gell-Mann matrix [TS] , then compute the integrals both in the coordinate 
and momentum spaces to obtain the correlation functions ^jLjHip)^ therefore the QCD spectral 
densities qcdA nnd qcdA nt the quark level through the dispersion relation. 
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In Eq.(18), we retain the term (gjcr^jyg'i) comes from the Fierz re-arrangement of the {qiA to 
absorb the gluons emitted from other quark lines to form {qjPsG'^pt'^^a^^qi) so as to extract the 
mixed condensate {qgscGq). 

Once the analytical QCD spectral densities Pjj^j^ qcdA qcdA obtained, we 

can take the quark-hadron duality below the continuum thresholds so and introduce the weight 
function exp (—^) to obtain the following QCD sum rules: 
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^ ds [Ap]L3H,QCDis) - P°3L3H,QCDis)] exp , 

( 22 ) 
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where pI,^,qcd(^) = ^QCois), 

PjLJH,QCD{s) = p]j^jH,ois) + p]^jH,3is) + + Pit 5 ( 5 ) + PjLjHfii^) + P^LJH^sis) 

+ Pil,JJI.9(®) + P^LjHPoi^) > 

^jLjH,QCDis) = + ^LJH,3(s} + + ^Lm,5(s) + 

~^^LjH,9(^) + ^LjH,10(^) > (23) 

the explicit expressions of the QCD spectral densities Pj^j„j(s) and with i = 0, 3, 4, 

5, 6 , 8 , 9, 10 are shown in the appendix. Here we introduce a negative sign in the definition 

PjLjH,QCD(^) = ^ 


^‘^P'jhiH QCoi^) to warrant positive spectral densities qcd('®)’ 

/^’^jLiH.QCois) > 0 . 


(24) 


In this article, we carry out the operator product expansion to the vacuum condensates up to 
dimension- 10 , and assume vacuum saturation for the high dimension vacuum condensates. 

We differentiate Eqs. (21-22) with respect to then eliminate the pole residues and 

obtain the QCD sum rules for the masses of the pentaquark states, 
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4 Numerical results 

We take the vacuum condensates to be the standard values {qq) = —(0.24±0.01 GeV)^, {qgsCrGq) = 
m^qq), ttiq = (0.8±0.1) GeV^, = (0.33GeV)"^ at the energy scale p = I GeV [lfl[T8]. The 

quark condensates and mixed quark condensates evolve with the renormalization group equation. 
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In the article, we take the 


MS mass mcirric) = (1.275 ± 0.025) GeV from the Particle Data Group |19) . and take into account 
the energy-scale dependence of the MS mass from the renormalization group equation. 
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where t = log^, 60 = 61 = &2 = ^ = 213MeV, 296MeV 

and 339 MeV for the flavors n/ = 5, 4 and 3, respectively [TO] . 

In this article, we study the pentaquark configurations consist of a light-diquark, a charm 
diquark, a charm antiquark, and resort to the diquark-diquark-antiquark model to construct the 
currents to interpolate the hidden-charm pentaquark states. The hidden charm (or bottom) five- 
quark systems qqiq 2 QQ could be described by a double-well potential. In the five-quark system 
qqiq 2 QQ, the light quarks qi and q 2 combine together to form a light diquark in color 

antitriplet. 


<7i + <72 -i' , 
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the Q-quark serves as a static well potential, which binds the light diquark to form a heavy 

triquark T^^q^g in color triplet, 


+ ^ (29) 

while the Q-quark serves as another static well potential, which binds the light quark q to form a 
heavy diquark in color antitriplet. 


q + Q —t VgQ , (30) 

where the i, j and k are color indexes. Then the heavy diquark in color antitriplet combines 
the heavy triquark 'Tq^q^g in color triplet to form a pentaquark state in color singlet. 

Such a doubly-heavy pentaquark state is characterized by the effective heavy quark masses Mg 
(or constituent quark masses) and the virtuality V = \/Mp — (2Mq)2 (or bound energy not as 
robust), just like the doubly-heavy four-quark states [9l |20l [21] |22l [23] . The QCD sum rules have 
three typical energy scales T^, we take the energy scale, = 0{T^), and obtain 

energy scale formula, 


fi = ,JmI - (2Mc)2 , 


(31) 


to determine the energy scales of the QCD spectral densities. In previous work [5], we take the 
value Me = 1.8 GeV determined in the diquark-antidiquark type tetraquark states [idO], and 
obtain the values /r = 2.5 GeV and /i = 2.6 GeV for the hidden charm pentaquark states Pc(4380) 
and Pc(4450), respectively. The energy scale formula works well. 

In this article, we choose the Borel parameters and continuum threshold parameters sg to 
satisfy the four criteria: 

1 . Pole dominance at the phenomenological side; 

2 . Gonvergence of the operator product expansion; 

3. Appearance of the Borel platforms; 

4 . Satisfying the energy scale formula. 

It is difficult to satisfy the criteria 1 and 2 in the QCD sum rules for the multiquark states. 
In the QCD sum rules for the hidden charm (or bottom) tetraquark states (or pentaquark states), 
molecular states and molecule-like states, the integrals 

^ dspgcD{s)exp(^-^^ , (32) 


are sensitive to the heavy quark masses mg, where the pgcois) denotes the QCD spectral densities. 
Variations of the heavy quark masses lead to changes of integral ranges of the variable ds 

besides the QCD spectral densities pgcois), therefore changes of the Borel windows and predicted 
masses and pole residues. In calculations, we observe that small variations of the heavy quark 
masses mg can lead to rather large changes of the predictions [HI HOI IHl HU HSj , some constraints 
are needed to specialize the heavy quark masses mg. 

Now we write down the definition for the pole contributions and use a toy-model spectral 
density to illustrate how to enhance the pole contributions. 
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Figure 1: The pole contributions with variations of the rric in the toy-model, where the perpen¬ 
dicular line corresponds to the MS mass mc(mc) = 1.275 GeV. 


with fc = 0, 1, 2, 3, 4, 5. The simple spectral density pqcd{s) makes sense, as we can simplify the 
calculation by taking the rough approximations ~ Am\ and fh\ = ~ 4m^, 

see the QCD spectral densities in the appendix. For the hidden-charm tetraquark states, fc < 4; 
for the hidden-charm pentaquark states, fc < 5. 

In Fig.l, we plot the pole contribution with variations of the c-quark mass rric for the typical 
Borel parameter = 3.5 GeV^ and continuum threshold parameter sq = 25GeV^. From the 
figure, we can see that the pole contribution decreases monotonously with the increase of the rric 
and fc. The MS mass mc{mc) = 1.275 GeV at the energy scale p = rric cannot lead to pole 
contribution > 50% for the hidden-charm pentaquark states as kmax = 5. A smaller mc{p) (or a 
larger energy scale p), for example, mc{p) — 1.1 GeV, can lead to the pole contribution > 50%. 
However, we cannot choose large energy scales freely to enhance the pole contribution, as the quark 
condensates and mixed condensates increase slowly but monotonously with the increase of energy 
scale, which slows down the convergent speed in the operator product expansion. In this article, we 
resort to the energy scale formula p = Mp — (2Mc)^ with the value Me = 1.8 GeV determined 
in the tetraquark states [5] to determine the energy scales of the QCD spectral densities, which 
works well in the QCD sum rules for the pentaquark candidates Pc(4380) and Pc(4450) [5]. 

In previous work [3120], we observed that the pole contributions can be taken as large as (50 — 
70)% in the QCD sum rules for the diquark-antidiquark type tetraquark states qq'QQ {X, Y, Z), if 
the QCD spectral densities obey the energy scale formula p = ~ The operator 

product expansion converges more slowly in the QCD sum rules for the pentaquark states qqiq 2 QQ 
compared to that for the tetraquark states qq'QQ. In Ref.[S], we observe that if we take the energy 
scale formula to determine the QCD spectral densities, the pole contributions can reach (40 — 60)%. 
So in this article, we try to choose analogous pole contributions, (50 ± 10)%. 

For the tetraquark states qq'QQ [3120], the Borel platforms appear as the minimum values, 
and the platforms are very flat, but the Borel windows are small, T^ax ~ = 0.4 GeV^, where 

the max and min denote the maximum and minimum values, respectively. For the heavy, doubly- 
heavy and triply-heavy baryon states qq'Q, qQQ', QQ'Q" [13 El, the Borel platforms do not 
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appear as the minimum values, the predicted masses increase slowly with the increase of the Borel 
parameter, we determine the Borel windows by the criteria 1 and 2 , the platforms are not very 
flat. The pentaquark states are special baryon states, as they have one unit baryon number. In 
this article, we also choose small Borel windows = 0.4 GeV^, just like in the case of 

the tetraquark states [a ED], and obtain the platforms by requiring the uncertainties induced 
by the Borel parameters are about 1%. In Ref.[S], we observe that such a criterion can be satisfied 
for the hidden-charm pentaquark states. 

Now we search for the optimal Borel parameters and continuum threshold parameters so 
according to the four criteria. The resulting Borel parameters, continuum threshold parameters, 
pole contributions, contributions of the contributions of the vacuum condensates of dimension 9 
and dimension 10 are shown explicitly in Table 1. From the Table, we can see that the criteria 1 
and 2 of the QCD sum rules are satisfied. 

In calculations, we observe that 

^ t i , 

Mi Mpt, (35) 

from the QCD sum rules in Eqs.(25-26). We can rewrite Eq.(31) into the following form, 

M|, = 4M2 -p ^2 ^ (36) 

which indicates that 

/it Mpf, 

/it Mp t . (37) 

It is difficult to obtain the optimal energy scales /i and masses Mp, however, the optimal energy 
scales /i and masses Mp do exist, see Table 2. 

We take into account all uncertainties of the input parameters, and obtain the values of the 
masses and pole residues of the ^ hidden-charm pentaquark states, which are shown in Figs.2- 
3 and Table 2. In Fig.2, we plot the masses with variations of the Borel parameters at large 
ranges, not just in the Borel windows. In the Borel windows, the uncertainties induced 

by the Borel parameters < 1%. From Table 2, we can see that the predicted masses have the 
relations Mqo.- < Moo,+ and Moi,_ < Moi,+ , which is consistent with our naive expectation, the 
pentaquark state with an additional P-wave has larger mass than corresponding S-wave state. The 
value Moi,_ = 4.30 ± 0.13 GeV is smaller than the value Mp^( 43 go) = 4.38 ± 0.13 GeV [5], which is 
also consistent with our naive expectation that additional unit spin can lead to larger mass. 

In the conventional QCD sum rules for the mesons, we usually take the continuum threshold 
parameters yso = Mgr -I- (0.4 — 0.6) GeV based on the assumption that the energy gap between 
the ground states and the first radial excited states is about 0.5 GeV, where the gr denotes the 
ground states. In Refs. fTBl E4] . we separate the contributions of the negative parity baryon states 
from that of the positive parity baryon states unambiguously, study the ^ and | heavy, 

doubly-heavy and triply-heavy baryon states qq'Q, qQQ', QQ'Q" with the QCD sum rules in a 
systematic way, the continuum threshold parameters yCo = Mgr -I- (0.6 — 0.8) GeV work well, the 
experimental values of the masses can be well reproduced. 

The pentaquark states are special baryon states, as they have one unit baryon number. In 
Ref. [5], we take the continuum threshold parameters ^so = Mpaiis&o/iibo) + (0.6 — 0.8) GeV, 
which also work well. In this article, the optimal continuum threshold parameters are = 
Mp -I- (0.6 — 0.8) GeV. One may worry that there maybe exist some contaminations from the high 
resonances and continuum states, as the spectroscopy of the pentaquark states is unclear in the 
present time. We should not be so pessimistic as the high resonances and continuum states are 
greatly suppressed by the factor exp (—^). If we take the largest threshold parameters and 
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Figure 2: The masses of the pentaquark states with variations of the Borel parameters T^, where 
the A, B, C and D denote the pentaquark states Too,-, Poi,-, Poo,+ and Poi,+, respectively. 
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Figure 3: The pole residues of the pentaquark states with variations of the Borel parameters 
where the A, B, C and D denote the pentaquark states Poo,-, Poi,-, Poo,+ and Poi.+ , respectively. 
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the central values of other parameters, then 



< 10 %, 


(38) 


the contaminations are greatly suppressed compared to the ground states, so the predictive ability 
cannot be impaired remarkably. The present predictions can be confronted with the experimental 
data in the future. 

In Fig.4, we plot the contributions of the pole terms with variations of the continuum threshold 
parameters and Borel parameters for the pentaquark states Poo,-, Poi,-, Poo,+ and Fbi,+ 
at the energy scales presented in Table 2. From the figure, we can see that the pole contributions 
decrease quickly and monotonously with the increase of the Borel parameters for the pentaquark 
states Poo,-, Poi,- and Fbo,+ ) the pole contributions reach 50% at « 3.3GeV^ with the central 
values of the continuum threshold parameters. For the pentaquark state Poi,+, the integral 

^ ds [yspoi, qcd(s) - Wc^i,qcd(s)] exp (-^) < 0. (39) 


at the value < 2.0 GeV^, which magnifies itself by the strange behavior of the pole contribution 
in Fig.4-D; while at the value > 2.3 GeV^, the integral is positive, the pole contribution decreases 
quickly and monotonously with the increase of the Borel parameter, and reaches 50% at the 
« 3.6 GeV^. We can draw the conclusion tentatively that the convergent behavior of the Poi,+ 
differs from that of the Poo,-, Poi,- and Poo,+ significantly, as it has much larger pole contribution 
in the Borel window, see Table 1. On the other hand, if we try to obtain smaller pole contribution, 
say about (40 — 60)% by choosing larger Borel parameters, the energy scale formula in Eq.(31) 
cannot be satisfied. From Fig.2-D, we can see that the Borel platform of the predicted mass 
Moi,+ appears as the minimum value, and the platform is very flat, which originates from the 
special convergent behavior in the operator product expansion. The negative integral at the value 
< 2.0 GeV^ or < 1.4 GeV shown in Eq.(39) is acceptable, as the optimal energy scale 
^ = 3.2 GeV 1.4GeV (see Table 2 or Fig.5-D), the value < 2.0 GeV^ or < 1.4 GeV is 
out of the allowed region of the Borel parameter = (3.0 — 3.4) GeV^, where the four criteria 
of the QCD sum rules can be satisfied. If we take into account the higher excited states besides 
the ground state, a larger continuum threshold sq is needed, therefore larger Borel parameter 
is needed to magnify the contributions of the higher excited states, then integral in Eq.(39) is also 
positive. So in the allowed region of the Borel parameter, the integral in Eq.(39) is positive. The 
continuum contributions can be approximated as 


^ dspH{s)exp(-^'^ 


where the ph{s) denotes the hadronic spectral density. At the value < 2.0 GeV^, the continuum 
contributions are greatly depressed, for example, exp (—fS-) < exp = 4.7 x 10“”^, and it 

is out of the allowed region of the Borel parameter. Furthermore, in the limit —>■ oo or in the 

local limit. 



ds ['/spoi^Qcoi^) 


'mc^oi,QCD{s)] exp 



-)> 

> 



ds \y/spoi,QCD{^) 


™c^1,QCd('S). 


0 , 


(41) 


a positive spectral density can be warranted. 
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T'(GeV') 



T'{GeV') 


Figure 4: The pole contributions of the pentaquark states with variations of the Borel parame¬ 
ters T^, where the A, B, C and D denote the pentaquark states Poo,-, Poi,-, Poo,+ and Poi,-i-, 
respectively. 



T2(GeV^) 

^/5^(GeV) 

pole 

(qq)^ 

(qgsCrGq)'^ 

-Poo,- 

3.1-3.5 

5.0 ±0.1 

(43 - 64)% 

(12 - 17)% 

(2 - 3)% 

Poi,- 

3.1-3.5 

5.0 ±0.1 

(40 - 62)% 

(13 - 18)% 

1 

CO 

Poo,+ 

3.1-3.5 

5.1 ±0.1 

(39 - 62)% 

(14-20)% 

1 

CO 

-Poi,+ 

3.0-3.4 

5.5 ±0.1 

(56 - 76)% 

-(6-12)% 

(3 - 6)% 


Table 1: The Borel parameters, continuum threshold parameters, pole contributions, contributions 
of the vacuum condensates of dimension 9 {{qq)^) and dimension 10 {{qggCrGq)'^). 



T2(GeV^) 

\/^(GeV) 

^(GeV) 

Mp(GeV) 

Xp{GeV^) 

Poo.- 

3.1-3.5 

5.0±0.1 

2.3 

4.29 ±0.13 

(1.39 ±0.26) X 10"^ 

P’01,- 

3.1-3.5 

5.0±0.1 

2.4 

4.30 ±0.13 

(2.36 ±0.45) X 10"^ 

Poo,-i- 

3.1-3.5 

5.1±0.1 

2.5 

4.41 ±0.13 

(0.60 ±0.12) X 10"^ 

Poi,-l- 

3.0-3.4 

5.5±0.1 

3.2 

4.82 ±0.08 

(3.11 ±0.37) X 10"^ 


Table 2: The Borel parameters, continuum threshold parameters, energy scales, masses and pole 
residues of the pentaquark states. 
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T'(GeV') 


Figure 5: The masses of the pentaquark states with variations of the Borel parameters and 
energy scales /i, where the A, B, C and D denote the pentaquark states Poo,-, Poi,-, Tbo,+ and 
Poi,+, respectively. 
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In this article, the contributions of the vacuum condensates of dimension-* with 

* = 0, 3, 4, 5, 6, 8, 9, 10 are defined by 


D 


3 LiH,i,+ 


D, 


/ 4 mg ds exp (-^) 

nil i'l^PjLm^QCDis) - mc^^jn,QCD(s)] exp (-^) 
iZl ds exp (-^) 

nil ds [\/sp]ij„,QCZ5(s) + ^o^jljh,QCd(Z exp (-#) 


(42) 

(43) 


which do not warrant the contributions Doo^t^-, Uqi,*,-, Doo^i^+ and Zloi,i,-i- have the same positive 
(or negative) sign, see Table 1. On the other hand, if we define the contributions of the 

vacuum condensates of dimension-* by 


D 


nZl W~^p]uhM - exp (-^) 


3 L 3 H,i,+ 


D 


SZl ds [V~SPljn,QCDis) - ^-f^3L3H,QCD{s)\ exp (-^) 

nZl ds [V^p]r.3H,ii^) + ^o^l3h4Z exp (-^) 




fZids \/spLh.qcd(s)+»^c^^,h,QCd(s) exp(-^) 


(44) 

(45) 


contributions of the terms '/sp]^jH *(«) are greatly enhanced compared to the terms *(s), 

which maybe lead to the contributions Hoo,*,-! _Doo,i,-i- and Zloi,i,-i- have the same positive 

(or negative) sign. However, we have to take into account the contributions of the high resonances 
and continuum states at the phenomenological side in the QCD sum rules. 

The correlation functions ip) can be written as 


^JlJh (p) 


^C„(p^/*)(C>„(/i)) 


I 


imlifj.) 


d^PQCoM 
s — 


r«o 




ds 


Pqcd{s,p) 


s — 


(46) 


at the QCD side, where the Cn{p^,p) are the Wilson coefficients and the {On{p)) are the vacuum 
condensates of dimension-**. At the energy scale /* ^ Aqcd, the short-distance contributions 
at > /*^ are included in the coefficients C„(p^,p), the long-distance contributions at 
are absorbed into the vacuum condensates {On{p))- The correlation functions ^jLjgip) are scale 
independent, 


which does not warrant 


= 


d r° j^^ Pqcd{s,p) 
dpJimlif,) S-p"^ 


—>■ 0 , 


(47) 


(48) 


due to the following two reasons in the present QCD sum rules: 

1 . Perturbative corrections are not available, the higher dimensional vacuum condensates are fac¬ 
torized into lower dimensional ones therefore the energy scale dependence of the higher dimensional 
vacuum condensates is modihed; 

2 . Truncations so set in, the correlation between the threshold 4*7*^(p) and continuum threshold 
So is unknown, the quark-hadron duality is an assumption. 

We cannot obtain energy scale independent QCD sum rules even if perturbative corrections 
are available, for example, in the case of the conventional heavy-light mesons |25j . but we have 
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typical energy scales which characterize the live-quark systems uudcc according to Eqs. (28-31) and 
serve as the optimal energy scales of the QCD spectral densities. In Fig.5, we plot the predicted 
masses with variations of the Borel parameters and energy scales n for the pentaquark states 
Poo,-, ^ 01 ,-, fbo,-(- and Poi,-i-, respectively. From the figure, we can see that the predicted masses 
decrease monotonously with increase of the energy scales /x. If we take the central values of the Borel 
parameters presented in Table 2, the uncertainties induced by the uncertainties 6fj, = ±0.6 GeV are 
about lo o 4 GeV, io;o 4 GeV, Ig^oa GeV and lo;o 4 GeV for the pentaquark states Poo.-, ^bi.-, -P 00 .+ 
and Poi,+, respectively. We can draw the conclusion tentatively that the uncertainties induced by 
the uncertainties of the energy scales in the vicinity of the optimal values are small. In calculations, 
we search for the optimal Borel parameters and continuum threshold parameters sq to reproduce 
the masses of the pentaquark states to satisfy the energy scale formula in Eq.(31). In other words, 
we take the energy scale formula in Eq.(3I) as a constraint, and do not take the energy scales of 
the QGD spectral densities as input parameters. 

The diquark-diquark-antiquark type current with special quantum numbers couples potentially 
to special pentaquark states. The current can be re-arranged both in the color and Dirac-spinor 
spaces, and changed to a current as a special superposition of the color singlet baryon-meson 
type currents. The baryon-meson type currents couple potentially to the baryon-meson pairs. 
The diquark-diquark-antiquark type pentaquark state can be taken as a special superposition of a 
series of baryon-meson pairs, and embodies the net effects. The decays to its components (baryon- 
meson pairs) are Okubo-Zweig-Iizuka super-allowed, but the re-arrangements in the color-space 
are non-trivial |26j . 

In the following, we perform Fierz re-arrangement to the currents Jqo ^^nd Jqi both in the color 
and Dirac-spinor spaces to obtain the results, 

Joo = cc — -Sc cu ± C 7 aC — —S^°^ccrfaU — -5(Tq,/ 3U ccr“^c ± -Saapcca°‘^u 

4 4 4 4 8 8 

1 1 It 

±-57“75MC7a75C- -57“75CC7a75U - -575MCZ75C± -Sj5CCi^5U , (49) 


Jqi 


-Succ- Sccu- -S-y^^ucjaC- -S'y°^ccjaU+ 2‘57“75MC7a75C± -57“75CC7„75W 
—iSj^u CZ75C — iSj5Ccij5U , (50) 


where we use the notations iSTc = C^^djTck and STu = C^^djTuk for simplicity, 

here the T denotes the Dirac matrixes. 

The components 5(a;)rc(a;)c(a;)r'u(a;) and 5(a:)rM(a;)c(x)r'c(a;) couple potentially to the baryon- 
meson pairs. The revelent thresholds are = 3.922 GeV, Mj/^p = 4.035 GeV, = 

4.151 GeV, = 4.293 GeV, M^^„p = 4.353 GeV, n{i 44 o) = 4.414 GeV, = 4.449 GeV, 

^A+( 2595 )bo = 4.457 GeV, = 4.463 GeV, M^+( 2595 )b.o = 4.599 GeV, = 4.604 GeV, 

D°{2420) ~ 4.708 GeV, -^a+B°(2430) “ 4.713 GeV [H]. After taking into account the currents- 
hadrons duality, we obtain the Okubo-Zweig-Iizuka super-allowed decays. 


Poo.-(4290) 
Poo.+ (4410) 
f’oi.-(4300) 
Poi,+ (4820) 


^ pJ/x/), A+P°, (51) 

^ pJ/V', A+P)*°,pryc, A+Z)°,pxc0, (52) 

^ pJ/xA, A+P*°,pry,, A+PO, (53) 

^ pJ/^P , A+P*° , A+(2595)P*° , ppe, V(1440)7, , pxci , A+pO , A+(2595)P° , 
A+P°(2420/2430), pxco , A+P°(2400), (54) 


where we add the masses of the pentaquark states in the brackets. We can search for the 
Poo.-(4290), Poo.+ (4410), Poi,_(4300) and Poi,+ (4820) in the those decays in the future. 


15 


5 Summary and discussions 

In this article, we present the scalar-diquark-scalar-diquark-antiquark type and scalar-diquark- 
axialvector-diquark-antiquark type pentaquark configurations in the diquark model firstly, then 
construct both the scalar-diquark-scalar-diquark-antiquark type and scalar-diquark-axialvector- 
diquark-antiquark type interpolating currents, and study the masses and pole residues of the 
— i hidden-charm pentaquark states in details with the QCD sum rules by calculating the 
contributions of the vacuum condensates up to dimension-10 in the operator product expansion. 
In calculations, we use the formula /i = Mp — (2Mc)^ to determine the energy scales of the QCD 
spectral densities. We can search for the pentaquark states Poo,-(4290), Poo,-i-(4410), Poi,-(4300) 
and Poi,-i-(4820) in the decays listed in Eqs.(51-54), and confront the present predictions of the 
masses to the experimental data in the future. 

The LHCb collaboration studied the A° —>• J/tppK~ decays and observed two pentaquark 
candidates Pc(4380) and Pc(4450) in the J/ipP mass spectrum [T]. The A° can be well interpolated 
by the current J{x) = {x)C"f^dj{x)bk{x) [24], the u and d quark in the A° form a scalar 

diquark [uc ?]3 in color antitriplet, the decays A° —>■ J/’tljpK~ take place through the mechanism 
A°([Md]6) —)> [ud\ccs —)> [ud]ccuus —>■ {[ud\[uc]c)K~(us) —>■ J/tf}pK~ at the quark level. We 
can also search for the pentaquark states Poo,-(4290), Poo,-i-(4410), Poi,-(4300) and Poi,-i-(4820) 
predicted in the present work in the decays A° —>■ J/^pK~ as the same mechanism works. 
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Appendix 

The QCD spectral densities and (with i = 0, 3, 4, 5, 6, 8, 9, 10) of the pen¬ 

taquark states. 


^ 4915207 r 8 / - V - zf [s - rnlf ( 8 s - 2,1711) , 

^o,o(s) = 9830407 r 8 / (v + z){l - y - zf (s - (7s - 2ml) , (55) 

Poi.o(s) = ^ 22807 r 8 / -y-zf{s- mlf ( 8 s - 3m2) , 

^i.o(s) = j dydz{y + z){l-y-{ls-2rril) , (56) 

Poo, 3 ( 5 ) = - ^ 72^6 J iy + z)ii-y- zf (s - mlf , 

^o. 3 ('S) = - ^ 36^6 / ^y^^ (1 - y - zf (s -mlf , (57) 

Poi. 3 ('S) = J dydz{y + z){l-y-zf (s-ml)^ , 

^i.3('S) = - J ^y^^ (1 - y - zf (s -ml)^ , (58) 
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^0,4('5) = 

Poi,4('^) = 

^ 1 , 4 ( 5 ) = 


TOt 


1474567r6 

19 ,asGG 


7864327r6 

13 .a,GG 


1310727r6 
TOg ,a 
' 1474567r6 
1 

'2949127r6 
19 


11796487r6 

13 MsGG 


J ^ (1 - y - (s - ™c) (5s - sto^) 

J {y + z){l-y- zf (s - (2s - ml) 

J dydzyz{l-y-zf (s-ml)^ {2s-ml) , 

(^ + ^ + ^ + ^) (1 - y - - ^c) ( 2 s - m2) 

J dydz (1 - y - (s - m2)" (5s - 2m2) 

^ UsGG ^ J (^i _ y _ 2)3 _ j7j2^2 _ 2rnl) 


(59) 


m“ ,a 


368647r6 

13 ,asGG 


327687r6 

m? .UgGG 


- 7864327 r 6 ^~^~^ J (?/ + ^((l - y - zf {s - ml)^ (5s - 2 m 2 ) , 

J (1 _ y _ 2)4 (s _ ^ 2 ^ (5s _ 3^2^ 

^ asGG ) J ^j^^ 2 yz(l-j/-z )2 (s - m^)" (2s - rn^) , 

+ ^ +J + (l-y-^)" (s-m 2 ) ( 2 s-rn 2 ) 

J (1 _ y _ 2)4 (s _ ^ 2^2 _ 2 — 2 ^ 

^ UgGG ^ J ^l_ y _ 2)3 _ 77^2^2 _ 2 _ 2 ^ 

^ asGG ) J _y_zf(^s- m2)" (5s - 2m2) , (60) 


737287r6 

1 


1474567r6 

19 ,a,GG 


5898247r6 

13 ,asGG 


3932167r6 
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Poo,5 (■5) — 


P00,5(®) — 


Poi,5('^) ~ 


= 


l9mc{qgsCrGq) 

327687r6 

mcjqgstjGq) 

40967r6 

l9mc{qgsCrGq) 

163847r6 

mcjqgsO-Gq) 

40967r6 

l9mc{qgscyGq) 
163847r6 
mcjqgsO'Gq) 
327687r6 ^ 

mcjqgsO'Gq) 
^ 491527r6 

ImcjqgsoGq) 

40967r6 J 

mcjqgsoGq) 
163847r6 ^ 

mcjqgsoGq) 
245767r6 


I dydz {y + z){l - y - z) {s - 
^ dydz (1 “ y “ (« -mlf , 

I dydz {1-y- z){s- ml)^ 

^ dydz (^ + ^{l-y-zf [s-mlY , (61) 

I dydz {y + z){l - y - z) [s - ml)^ 

^ dydz + 

^ dydz + (® “ ^c) (5s “ 3mc) , 

dydz {1 — y — z) [s — 

^ ^y'^^ (y z) ~ ^ ~ ~ 

^ dydz (1 — y — z)^ (s — mf) (2s — mf) , (62) 


Poo,6(s) 

^o,6(s) 

Poi,6(s) 

(^l,6(s) 


I ~ y ~ {s — oic) (5s — 3TOg), 

J dydz (y + z)(l - y - z) {s - ml) (2s - ml) , 


1927r4 

Ml 

1927r4 

M! 

487r4 

M! 

967r^ 


j dydzyz(l-y-z) (s-ml) (5s - 3^^) , 

J dydz jy + z)(l - y - z) {s - ml) (2s - ml) , 


(63) 


(64) 


1 f N ^Hqq){qgsoGq) f , 2\ 

Poo,sis) =-- dydzyz(4s-3m^) 


P00,8(s) — 

Poi,8(s) = 
Poi,8(s) = 


{qq){qgsoGq) 

40967r4 

(3, 

{qq){qgsoGq) j 

20487r4 J 

55i{qq){qgsoGq) 

15367r4 

(3,-2M?) 

jqq)jqgsoGq) 


dydz jy + z)(l - y - z) (4s - iml) , 


dydz j\ — y — z) (3s — 2ml^ , 
f dydz yz (4s — Sfnl) , 


10247r4 


dydz jl — y — z) (3s — 2ml) > 


(65) 


( 66 ) 
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Pooao('®) = 

^o,io('®) = 

Poi,io('^) = 


Poo ,9 (' 5 ) = 

mdqq)^ ^ 

- 144 ,^ i 


^ 0 , 9 (' 5 ) = 

mdqqf ^ 

- 1 , 

( 67 ) 

Poi, 9 ('®) “ 

mdqq)^ ^ 

- i 


Po 1 , 9 (' 5 ) ~ 

mdqq)^ ^ 

- 18 .^ 1 , 

( 68 ) 

19 (ggsCrGg)^ 

245767 r 4 J 

rvf 

1 dyy{l-y )[3 + ml 6 {s-ml)] 

Vi 


17 {qgsaGq)'^ 

1474567 r 4 

■ J dydz {y + z) [3 + <5 (s - m^)] , 


19 {qgsaGqy 
491527 r 4 j 

/•Vf 

1 dy [2 mis (s — ml)] 

Vi 


17 {qgsaGq)'^ 

737287 r 4 

J dydz [2 + (5 (s — m^)] , 

( 69 ) 

19 (qgsCrGq)'^ 
61447 r 4 j 

rVf 

1 dyy{l -y )[3 + mlS{s- ml)] , 

Vi 


19 (qgsaGq)'^ 
245767 r 4 j 

fVf 

1 dy [2 + mlS [s- ml)] 

Vi 


17 {qgsaGq)'^ 

368647 r 4 

J dydz [ 2 -\-rnl 5 [s— ml)] , 

( 70 ) 


y^l ^2 _ (?/+z)™c 


Zi = — — m~ 

'' ys—m^ ‘ ^ yz 


where / dydz = dy dz, yf = y* = ^ 

= y(i-y) ’ ly^ fo fz- ^ Jo ^ when the 5 functions <5 (s — m^) and 5 (s — m^) 
appear. 
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